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ABSTRACT 

We present mid-infrared Spitzer IRAC and MIPS images of the UX Orionis star VV Ser and the 
surrounding cloud. The 5.6-70 /itm images show bright, localized and nebulous emission extended 
over 4' centered on VV Ser. We interpret the nebulosity as being due to transiently heated grains 
excited by UV photons emitted by VV Ser. A companion paper describes how the physical structure 
of the VV Ser disk has been constrained using a wide range of observational data modeled by an 
axisymmetric Monte Carlo radiative transfer code. The model of the system in particular constrains 
the strength of UV radiation field throughout the envelope, and the nebulosity thus presents an 
opportunity for directly constraining the excitation mechanism and opacity law of PAHs in molecular 
clouds. In this paper we employ the model to study the nebulosity surrounding VV Ser using quantum- 
heated PAH molecules and Very Small Grains (VSGs) consisting of amorphous carbon in the thermal 
cooling approximation. Imprinted on the nebulosity is a wedge-shaped dark band, centered on the 
star. We interpret this dark wedge as the shadow cast by the inner regions of a near-edge-on disk 
in UV light, allowing the PAHs to be excited only outside of this shadow. The presence of a disk 
shadow strongly constrains the inclination as well as the position angle of the disk. The nebulosity 
at 5.6-8.0/xm and the 2175 A absorption feature seen in an archival spectrum from the International 
Ultraviolet Explorer can be fit using only PAHs, consistent with the main carrier of the 2175 A feature 
being due to the graphite-like structure of the PAHs. The PAH component is found to be relatively 
smoothly distributed in the cloud, while the population of VSGs emitting at 20-70 /xm is strongly 
concentrated ~ 50" to the south east of VV Ser. Although depending on the adopted PAH opacity, 
the abundance of PAHs in the surrounding cloud is constrained to 5 ± 2% of the total dust mass, 
given the opacity. The extent of the nebulosity constrains the density of the gas surrounding the VV 
Ser disk to 500 ± 200 cm~ 3 for a gas-to-dust ratio of 100. This low density suggests that the quantum 
heated material is not part of the original envelope of VV Ser and that it is rather a quiescent part 
of the Serpens molecular cloud that the star has passed into after being formed. Although relatively 
rare, quantum heated nebulosities surrounding single, well defined stars are therefore well-suited for 
gaining unique insights into the physics of very small particles in molecular clouds. 
Subject headings: accretion, accretion disks - circumstellar matter - stars: formation, pre-main- 
sequence - infrared: stars 



1. INTRODUCTION 

It is well known that if a sufficiently small grain is 
hit by a high energy (ultraviolet) photon, the internal 
energy cont ent far exceeds that of a grain in thermal 
equilibrium ijDulevll 19731 ISe llgrcn 1984). Once the small 
grain has been heated, it will rapidly cool by emitting 
strongly in the mid-infrared wavelength regime. Yet, 
because the energy is distributed in many vibrational 
modes (typically 3(V atom — 2)), the subsequent cooling 
of the grain takes place in an almost thermal (classi- 
cal) manner. Such quantum-heated grains have been 
used to explain a range of observations in the past. 
Specifically, emission from polycyclic aromatic hydrocar- 
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bons (PAHs) require quantum-heati ng l|Leger fc Pugetl 
119841 ISiebenmorgen fc K rucgcl 1992), but also the pres- 
ence of very small carbon/silicate grains (VSGs) of sizes 
5 — 50 A has been inferred by infrared o bservations 
ijBoulanger et all 119881: iHenning et aTlll998f) . In many 
extragalactic environments, observations with the In- 
frared Space Observatory (ISO) have shown the impor- 
tance of emission from quantum-heated grains. VSGs 
have been suggested to play a crucial role for the emission 
spectra of especially starburst galaxies llContursi et alJ 
|200H iForster Schreiber et all 12001 iRouan et al.T T2004 h 
Typical for all environments with high UV fields is that 
PAHs dominate the spectrum from 5-15 /mi, while VSGs 
dominate at wavelengths > 15 //m. 

This paper is the second part of a study of Spitzer 
observations of the Herbig Ae star VV Ser. We fo- 
cus on VV Ser for several reasons. Mid-infrared im- 
ages obtained with the Spitzer Space Telescope as part 
of the Legacy Pro gram "From Molec ular Cores to Proto- 
planetary Disks" l|Evans et al.ll2003(l show a bright bipo- 
lar nebulosity surrounding the central star. This was a 
considerable surprise since the nebulosity has no or lit- 
tle optical and near-infrared counterparts and therefore 
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the density must be low such as not to create any de- 
tectable reflection nebulosity. The simplest way of ex- 
plaining this, is if the nebulosity observed by Spitzer is 
caused by very small (< 0.01 /im) transiently heated dust 
grains. 

Another important reason for selecting VV Ser for fur- 
ther study is the presence of a dark band bisecting the 
entire mid-infrared nebulosity. The dark band has a 
"wedge-like" shape suggesting the presence of a nearly 
edge-on circumstellar disk. Yet, the extent of the band of 
> 4' or 60,000 AU is far too large for a circumstellar disk. 
Exactly this type of morphology has been seen in reflec- 
tion nebulosities sur r ounding a number of young star s 
ijHodapp et alJ 120041 iPontoppidan fc Dullemondl 12005^1 . 
These studies suggest that if an edge-on or nearly edge- 
on disk is surrounded by large-scale envelope material 
(or diffuse material), then the shadow cast by the disk 
into/onto the surrounding material in scattered light can 
be seen as a dark wedge-like band originating from the 
star and extending out to tens of thousands of AU, i.e. 
many times the size of the disk. There are good reasons 
to believe that any extended PAH emission excited by 
a single source should exhibit a similar behaviour due 
to the shadowing of ultraviolet photons by a circumstel- 
lar disk. In fact, the effect may be more pronounced 
in the case of quantum-heated grains since the exciting 
UV radiation from the central star will be subject to 
a stronger shadowing effect than near-infrared photons. 
The opening angle of a disk shadow and the sharpness 
of the shadow may constrain the structure and compo- 
sition of the disk as well as the opacity spectrum of the 
quantum-heated grains. Disk shadows are therefore a 
powerful tool for studying the physical structure of spa- 
tially unresolved circumstellar disks. Clearly, the best 
way of analyzing these data is the simultaneous fitting of 
all available data with radiative transfer models of disks 
surrounded by an envelope. 

There is considerable additional evidence that the 
VV Ser system is highly inclined (> 70°). This in- 
cludes broad, double-peaked CO rovibrational emis- 
sion lines from the fu ndamental band around 4.7/im 
ijBlake fc Booeertl 12004) . Near-infrared interferometric 
measurements, although som ewhat ambiguous, s eem t o 
point to a high inclination i|Eisner et all 1200.1 (2004). 
Most importantly, perhaps, VV Ser is a well-known UX 
Orionis variable star, i.e. the star experiences frequent, 
but irregular, brightness dips, presumably due to orbit- 
ing dust clumps in the disk p assing in front of the star . 
The first part of this study l|Pontoppidan et all 120061 
henceforth Paper I) describes the observations of VV Ser 
in terms of its status as a UX Orionis type variable and 
argues extensively for a high inclination of the VV Ser 
disk. In Paper I we construct an axisymmetric Monte 
Carlo radiative transfer model of the VV Ser disk and 
use the model to show that the Spectral Energy Distri- 
bution (SEP), the near-inf rared interferometric visibil- 
ities of lEisner et al.l il2()() 1 ■ as well as the UX Orionis 
variability can be explained in terms of a self-shadowed 
nearly edge-on disk. In the present study we use the 
same model with the addition of quantum heated grains 
to study the large extended nebulosity surrounding VV 
Ser. Our goal is to construct a single radiative trans- 
fer model that explains all the observed characteristics 
of VV Ser in terms of an inclined disk embedded in a 



relatively low density cloud containing PAHs and VSGs. 
Once the structure of the system as well as the radiation 
field is known, the physical properties of the quantum 
heated grains surrounding VV Ser can be constrained, 
in particular their abundance, and to some degree, their 
opacities. 

This article is organized as follows: In Sect. [3 we de- 
scribe the Spitzer imaging observations of VV Ser and in 
Sect. El the cloud environment in which VV Ser is em- 
bedded. Sect, ^describes the radiative transfer model 
and our treatment of quantum-heated grains. Sect. 
discusses how the Spitzer images and other imaging data 
compare to the model and Sect. El provides a discussion 
in particular of how quantum-heated nebulosities around 
isolated low-mass stars can be used to constrain the small 
grain population associated with star formation. 

2. OBSERVATIONS 

Most of the observations of VV Ser used in the mod- 
eling are described in Paper I. The mid-infrared spec- 
troscopy and ima ging were obtained with the Spitzer 
Space Telescope l)Werner et al.l l2004|) as part of the 
Legacy program, "F rom Molecular Cor es to Protoplan- 
etary Disks" (c2d) (lEvans et all 120031) using all the 
available in s trume nts [IRS, iHouck et al.l 120041) . I RAC , 
iFazio et alJ <|2004[) and MIPS. iRieke et all l|200l ]. In 
addition to the Spitzer-IRS mid-infrared spectroscopy of 
VV Ser presented in Paper I we use IRAC and MIPS 
imaging at 3.6, 4.5, 5.6, 8.0, 24 and 70 /mi of the area 
surrounding VV Ser. The im aging was reduced u sing the 
c2d mosaicking pipeline (see lHarvev et alJ l)2006tk Spies- 
man et al. in prep, for IRAC and MIPS, respectively). 
In support of the Spitzer observations, we have also ob- 
tained a continuum map at 850 ^m with SCUBA on the 
James Clerk Maxwell Telescope (JCMT) (see also Paper 
I). From the European Southern Observatory archive, 
we have extracted archival images in the J-, H- and K- 
bands obtained with SOFI at the New Technology Tele- 
scope (NTT) 6 . The near-infrared images were reduced 
using standard methods of dark subtraction, flat field 
correction and registering of individual frames using cross 
correlation. An ultraviolet spectrum obtained with the 
International Ultraviolet Explorer (IUE) was obtained 
from the IUE archive and is described in Paper I. 

3. THE MOLECULAR SURROUNDINGS OF VV SER 

VV Ser was identified as a particularly interesting 
object when mid-infrared IRAC and MIPS images re- 
vealed that the star is surrounded by a bright, extremely 
large nebulosity not seen in near-infrared images. An 
IRAC/MIPS color composite image is shown in Fig. ^ 
The nebulosity extends over an area at least 6 arcmin- 
utes (or 94,000 AU at a distance of 260 pc) across and is 
brightest in the 8.0, 24 and 70 /jm bands (IRAC4, MIPS1 
and MIPS2), but is also weakly detected in the 5.6/mi 
band (IRAC3). At 3.6 and 4.5 /jm (IRAC1 and IRAC2), 
no extended emission is detected. We interpret the neb- 
ulosity as being due to a quantum-heated component of 
either an envelope or simply the surrounding molecu- 
lar cloud. The relatively low extinction (Ay ~ 3 mag) 
toward VV Ser strongly suggests that little original en- 
velope material is left. This quantum-heated component 

6 In part based on observations obtained at the European South- 
ern Observatory at La Silla under program 67.C-0042(A). 



-4xl0 4 -2xl0 4 2xl0 4 4xl0 4 -4xl0 4 -2xl0 4 2xl0 4 4xl0 4 

AU AU 

Fig. 1. — 4.5 (blue), 8.0 (green), 24.0 fi (red) color composites of the VV Ser nebulosity. Left panel: Spitzer-IRAC bands 2 and 4 together 
with Spitzer-MIPS band 1. The field is centered on VV Ser and measures 400" on the side. The bright red source to the left of VV 
Ser is a presumably unrelated young star, a background evolved star or possibly a very bright background infrared galaxy. In this color 
scheme, background stars are blue, while PAH emission is green. Right panel: Model color composite. The model has been convolved 
with theoretical PSF profiles calculated using the STINYTIM package ver. 1.3 for pre-launch parameters (Krist 2002, SSC web site). 
Additionally, Gaussian noise has been added to the model with a standard deviation corresponding to that of the actual Spitzer images. 



likely consists partly of Polycyclic Aromatic Hydrocar- 
bons (PAHs) and partly of Very Small Grains (VSGs). 
The presence of the VSGs is inferred from the bright 24 
and 70/Ltm nebulosity (see Sect. I5.1|) . The presence of 
such compact, mid- infrared nebulosities is clearly a rare, 
but not unique, phenomenon, judging from the larger c2d 
Spitzer maps. 

The optical colors of the star correspond to a steady 
extinction of Ay ~ 3, but with frequent, non-periodical 
dips lasting of order 10 days of 0.5-3.0 additional magni- 
tudes of extinction. Using the 2MASS point source cat- 
alog, an average extinction through the Serpens cloud in 
a 5' x 5' field around VV Ser is determined by fitting 
a reddened main sequence to the colors of background 
stars in the field. Assuming a near-infrared extinction 
law of A\ oc A -19 (Kaas 1999) an average extinction of 
Aj = 1.5 ± 0.2 is found. For an Ry = 3.1 optical ex- 
tinction law, this correspond s to an Ay ~ 6. T he lower 
resolution extinction map of lGa,mbresvi lfl999l) finds an 
extinction of Ay = 7 — 8 mag in the region. For a con- 
stant density molecular cloud, this places VV Ser roughly 
in the middle (along the line of sight) of a filament aligned 
in the north-south direction. The northern end of the fil- 
ament culminates in the famous Serpens Cloud Core [in- 
cidentally, the Serp ens Core contains CK 3, ano t her dis k 
shadow candidate, Pontoppidan & Dullemond (2005)1. 
The width of the filament is roughly 22' l)Cambresvl 
1999). Using the conversion bet ween column densit y 
and extinction in the J-band from IVuong et al.l l|2003j) , 
N(R 2 ) = 2.8 x 10 21 cm~ 2 x Aj, the molecular cloud den- 
sity around VV Ser is estimated to ~ 700 cm~ 3 . 

4. MODEL 

To model the observed SED and Spitzer imaging of 
VV Ser, we use the axis ymmetric Monte Carlo radiat ive 
transfer code RADMC ijDullemond fc Dominikll2004ft in 
combination with the rayt racing capabilities of the mor e 
general code RADICAL IjDullemond fc Turollal l2000j) . 
The density structure is axisymmetric, but the pho- 
tons are followed in all three dimensions. The RADMC 
code is used to derive the temperature structure and 



scattering source function of a given dust distribution. 
Once these parameters are known, images can be cal- 
culated using RADICAL. This set up has been used 
to model similar protostellar disks llPontoppidan et alJ 
l2005UPontoppidan fc Dullemondll2005|) . In Paper I it is 
described how the physical disk structure is constrained 
to fit the observed SED, images and UX Orionis be- 
haviour. In this paper, we add a simple treatment of 
quantum-heated grains to the model in order to describe 
the mid-infrared nebulosity surrounding VV Ser. We 
note that while the treatment of quantum heated grains 
included in the model is intended for the VV Ser neb- 
ulosity, it can easily be adapted to treat a wide range 
of axisymmetric radiative transfer problems. One obvi- 
ous possibility is the treatment of emission from PAH 
molecules in the disk itself, although that would require 
the inclusion of multi-photon excitations. The code only 
treats isotropic scattering, which is not a good approx- 
imation for the UV photons considered in the models. 
However, in terms the quantum heated nebulosity we ar- 
gue that the scattering phase function only plays a minor 
role. Since the optical depth to scattering of UV photons 
is only ~ 0.5 through the envelope, the majority of pho- 
tons heating small grains in the extended nebulosity have 
experienced at most a single scattering event from the 
disk surface. This introduces an error of up to a factor 
of two in surface brightness - similar to the uncertainty 
on the opacities used. 

4.1. Quantum heated grains and multi- dimensional 
radiative transfer 

Quantum heated grains are treated in the model 
using the "continuous coo ling" approximation (e.g. 
iGuhathakurta fe Draindll989|) . This popular approxi- 
mation is valid for sufficiently large grains/molecules. In 
principle, it should only be used for grains with more 
th an 100 atoms , but i n the case of PAHs it was found 
by IDraine fc Lil (|2001^ to provide accurate results for 
molecules as small as Csoil x by comparing to the full 
statistical treatment. The continuous thermal cooling 
approximation is a method for calculating PAH emission 
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spectra without having to treat the radiative transfer in 
a full quantum mechanical framework. When includ- 
ing quantum heating in a multi-dimensional radiative 
transfer code, this approximation results in an signifi- 
cant advantage in terms of CPU time relative to a code 
calculating a full transition matrix treating all upwards 
and downwards transitions between vibrational modes. 
We emphasize that the code presented here is also very 
portable to any kind of axisymmetric structure. 

In the adopted approximation, grains are heated in 
a quantized fashion, but cooled classically. The neces- 
sary assumption is that a sufficient number of vibra- 
tional modes are excited in the grain to create a con- 
tinuous cooling behaviour, at least in the temperature 
range where most of the power is radiated. In terms of 
the radiative transfer, a necessary assumption is that the 
absorbed and emitted radiation fields are decoupled. In 
the case of PAHs, this amounts to assuming that infrared 
photons do not excite the PAH molecules to significant 
temperatures. In the current version of the code, only 
single photon excitations are treated. Multi-photon exci- 
tations are important for very strong radiation fields and 
do not apply to the nebulosity surrounding VV Ser at dis- 
tances of 5,000-50,000 AU from the central star, where 
no additional photons will strike a given PAH molecule in 
the time it takes to cool to the background temperature. 
For inclusion of disk PAHs, multi-photon excitations will 
be essential and therefore we do not attempt to model 
the disk PAH features seen in the Spitzer spectrum. 

4.2. The thermal cooling approximation 

The treatment of quantum heated grains can in the 
thermal cooling approximation be regarded as a post- 
processing add-on to the radiative t ransfer Monte Carlo 
code. The Monte Carlo code uses theiBiorkma n fc Woodl 
i200H) method w ith continuous absorption according to 
the Ibucvl Ipfll algorithm. This method follows pho- 
ton packages from the surface of the central star as they 
travel through a gridded density structure. Photon pack- 
ages can be scattered, changing their directions, or ab- 
sorbed, altering both their wavelength and direction. 
The code supplies the UV radiation field at each grid 
point, and the calculation of the non-thermal emission is 
a matter of calculating a response function that can con- 
vert the fraction of UV photons absorbed by non-thermal 
grains into infrared photons. 

In essence, to avoid the computationally intensive task 
of re-calculating the temperature distribution function of 
the quantum heated grains every time a photon interacts 
with a grid cell, we split the problem in the following way: 
The mean intensity of quantum heating photons is ob- 
tained by adding the energy absorbed from each photon 
packet passing through a given cell. The mean intensity 
for each grid point is then saved and used at the end of 
the Monte Carlo run to calculate the source function of 
quantum heated grains, as described below. To facilitate 
this approach, the photon packages are divided into pri- 
mary and secondary components. The primary photon 
packages are those that originate at the stellar surface, 
and are only scattered. If a primatry photon package 
is absorbed by a quantum-heated grain, it is lost to the 
Monte Carlo code. If it is absorbed and re-emitted by 
a thermal grain, it is converted to a secondary photon 
package. All grains, including quantum-heated grains, 



are treated as thermal by a secondary photon package. 
After completing the Monte Carlo run, the primary pho- 
ton packages absorbed by the quantum-heated grains are 
used by an external code to calculate a quantum-heated 
emissivity. When this is done, a second Monte Carlo 
run can be executed in which photons can be launched 
from the quantum heated grains according to the emis- 
sivity calculated by the external code, in addition to pho- 
tons launched from the stellar surface. Photon packages 
launched by quantum heated grains begin as secondary 
grains, and will only heat thermal grains. Typically, the 
effects of this second Monte Carlo run on the SEDs and 
images are minor. 

For PAHs, this split between primary and secondary 
photon packages is likely accurate since the PAH opac- 
ity drops by at least 3 orders of magnitude above 1 /xm, 
according to the opacities of iLi fc Draind |2001). So 
even though there are 10-100 times more photons at 
near-infrared wavelengths than at UV wavelengths, only 
1-10% of the heating events will potentially be due to 
infrared photons. However, recent laboratory results 
have shown that at least some classes of ionized PAHs 
may have significant cross sections in the near-infrared 
ijMattioda et alJl2005j) . Clearly, future work will need to 
explore the effects of this in further detail. In terms of 
the model, an under-estimate of the near-infrared cross 
section may result in an over-estimate of the abundance 
of quantum heated grains. Note, however, that single in- 
frared photons will not be able to excite PAHs to the tem- 
peratures responsible for the short-wavelength PAH fea- 
tures below 10 /im. For carbonaceous or silicate VSGs, 
this assumption may be less valid because these types of 
grains have a less steep opacity curve. 

In the following, D is a frequency of the UV radiation 
field (primary) photon as opposed to ^, which indicates 
a frequency of a re-emitted secondary photon. 

The cooling rate of the quantum heated grains, as- 
sumed to be radiating thermally, is: 

poo 

r coo i= / ^B v {T)K v dv, (1) 

Jo 

where k„ is the opacity of the quantum heated grains, 
and B V (T) is the Planck function. The time to cool to a 
temperature, T, can then be calculated as: 

Tcooi(T) = / Tp ° ak ^^dT, (2) 

JT 1 cool U ) 

where C(T) = dU /dT is the heat capacity of the dust 
material. U(T) is the enthalpy curve related to the bulk 
material of which the quantum heated grains consist. 
The grain contains N^tom atoms, which in the case of 
a PAH molecule refers to the number of carbon atoms. 
The peak temperature, T poa k, reached when the grain 
is hit by a photon of energy v is simply found via the 
enthalpy curve: 

^(Tpcak) = hi>/N atom (3) 

The next step is to calculate the time-integrated emission 
per mass of the quantum-heated grains. This produces 
the output energy spectrum emitted if one photon of en- 
ergy hv is absorbed. 

Ev{V) = k v B v {T) x dT, (4) 
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The emissivity of the quantum heated grains is then de- 
termined by the UV radiation field alone, \v ( e -g- m 
units of ergs -1 cm -2 Hz -1 sterad -1 ), and the mass den- 
sity of transiently heated particles, pqh- Integrating over 
the photons absorbed by the quantum heated grains de- 
fines a response matrix, yielding an emissivity j v as a 
function of an input radiation field spectrum: 



Pi?KvpQB.-r^:E u {v)dv (5) 



We have included the effect of destroying the grains 
if Tp ea k is higher than the sublimation temperature. In 
essence this works as an energy sink since the absorbed 
energy is assumed to be used to heat the gas and subse- 
quently be emitted as gas-phase lines in an optically thin 
part of the spectrum. The grain destruction is described 
by the step function, [3 V : 



, = i 1 for T peak < 1500 K 
' v ' otherwise 



(6) 



In fact, for the single photon excitation case, no grains 
are destroyed, since a 0.1 /im photon will excite a C50H20 
PAH molecule to just under 1500 K and a 1000 atom 
silicate grain to 200 K. However, for the multi-photon 
excitation case, grain destruction may play an important 
role for the radiation fields present within a few AU of a 
Herbig AeBe star. 

4.3. PAH and VSG opacities and enthalpy curves 

The physics of a given population of quantum heated 
grains is contained in the opacity, in the enthalpy curves 
and the grain size. The presence of strong 24 and 70 /jm 
nebulosity in addition to the 6-8 /im emission suggests 
that both a PAH component (providing the short wave- 
length opacity) as well as a carbon and/or silicate VSG 
component (to provide the far-infrared opacity) are re- 
quired. The PAH opacity, kpah, has been calculated 
using the code of Visse r et al . (in prep), following the 
method of iDraine fe Lil 1)200 1|) . This opacity includes a 
graphite-like contin uum opacity and is similar to that of 
iLi fe Draind l)2001|) (see their Fig. 2). For this paper, we 
opted to use a single generic PAH molecule consisting of 
neutral Cso^o- In reality, the PAH component is much 
more complicated and requires the use of a full chem- 
ical code to determine. For the purposes of fitting the 
broad-band photometry of the nebulosity around VV Ser, 
a simple PAH opacity is deemed sufficient. The model 
framework developed for VV Ser is a good starting point 
for including a full PAH chemical network. 

The nature of the VSG component is even less con- 
strained. It is clear that a quantum-heated component 
is needed to produce the 24-70 ^m flux. To achieve 
this, the grains must be heated by UV-near-infrared pho- 
tons to temperatures of order 150-300 K. If UV photons 
are responsible, then grains significantly larger than the 
PAHs are needed. In the literature, both silicate and 
carbonaceous VSGs have been proposed, presumably as 
amorphous solids t o avoid producing strong spec tral res- 
onances. Following lWeingartner fe Draind l )2001|) . we as- 
sume that the VSGs are mostly carbonaceous in nature, 
with a size distribution corresponding to a small grain 
tail of the larger grains. For the VSGs, we use the opac- 
ities of spherical graphite grains in the Rayleigh limit. 
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Fig. 2. — The opacities used for the envelope. The PAH and 
VSG opacities has been scaled to 5% and 6% of the total dust mass, 
respectively. The solid red curve is the sum of the opacitie s . The 
total opacity is compared to that of Wcineartncr & Draine 1120011, 
for Ry = 5.5, which also takes PAH and VSGs into account. 



TABLE 1 

VV Ser envelope model parameters 



Parameter 



Model value 



Pout (5 x 10 4 AU) 500 ±200 cm" 3 
R in 15 000AU 
Rout 50 000AU 
PAH size C50H20 
VSG size 500 C 
Apah 5 ± 29c 1 
Avsg 6±3% 1 ' 2 
Disk shadow pos. angle 15 ± 5° 
Disk inch 72 ± 5° 

1 Abundances are given in terms of percentage of the total 
dust mass (quantum heated + thermalized grains). 

2 The VSG abundance is for the south-eastern clump. 

The enthalpy curve for the VSGs is the same as that of 
the PAHs with a single size of 500 C atoms. The exact 
choice in size makes little difference in the absorption 
coefficient, but is important when calculating the heat 
capacity for a single grain; smaller grains can be heated 
to higher temperatures by a photon of a given energy. 
Since the SED of the nebulosity is not well-constrained 
in the far-infrared, the size of the VSGs cannot be di- 
rectly constrained. For the PA H/VSG enthalpy c urve, 
we use the parametrized form of lChase et all II1985I) . rel- 
evant for graphite. The properties of the quantum heated 
grains are summarized in Tableland the opacity curves 
of the various envelope dust components are compared 
in Fig. □ 

5. FITTING THE ENVELOPE STRUCTURE 

The constructed model is able to fit data spanning four 
orders of magnitude in wavelength and five order of mag- 
nitude in spatial scale (see Paper I for a discussion of 
the fit on small spatial scales). In this section, it is dis- 
cussed how each observational property is connected to 
the model. 

5.1. Mid-infrared nebulosity 

A central property of the model of VV Ser is the 
ability to reproduce the bright nebulosity observed at 
mid-infrared wavelengths with Spitzer-IRAC and MIPS. 
Clearly, a mid-infrared nebulosity extending over 6' asso- 
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Fig. 3. — Spectral energy distribution of VV Ser in a 50,000 AU 
aperture (dashed green curve) and a 600 AU aperture (full red 
curve) the quantum-heated nebulosity. The □ symbols indicate 
IRAS fluxes from the IRAS extended source catalog, while the 
O symbols indicate Spitzer-MIPS fluxes for the entire nebulos- 
ity (upper) and the central VV Ser point source (lower). The 
JCMT-SCUBA point indicates an upper limit for a point source at 
the location of VV Ser. The model SED including the quantum- 
heated nebulosity has been calculated by integrating the flux over 
the south-eastern lobe within 50,000 AU, where most of the enve- 
lope emission is, and subtracting the background cloud emission 
measured in an annulus spanning 50,000-100,000 AU. 



dated with a young star is not a common occurrence. As- 
suming a distance of 260 pc mid-infrared emission is de- 
tected up to 47,000 AU from the central star. In IRAC4, 
MIPS1 and MIPS 2 (centered on 24 pm and 70 /im) a 
maximum surface brightness of 3, 4 and 18 mJy/sq. arc- 
sec is reached at a distance of 15,000AU (see Fig. EJ. 
This corresponds to color temperatures between IRAC4 
and MIPS1 of 360 K and between MIPS1 and MIPS 2 of 
85 K. Tests using the model showed that it is not possible 
for dust grains in thermal equilibrium with the radiation 
field from the central star to produce such bright nebu- 
losity in the mid-infrared at these radii, except at 70 fjxa 
where small thermal grains can explain the emission. The 
inclusion of quantum-heated grains readily produces a 
nebulosity with high surface brightness around a star 
as luminous as VV Ser. A possible way of producing 
non-thermal continuum emission is via quantum-heated 
VSGs or PAHs or a mix of both. In this context it is im- 
portant to note that the nebulosity is bright at all wave- 
lengths beyond ^5.6 /im (IRAC 3) and in fact continues 
to rise to dominate the total SED at 70 /im. It is also 
interesting to note that the 24 and 70 /im emission is not 
smoothly distributed around the source, but is strongly 
concentrated 50" to the south-east of the star. The 8/im 
emission, on the other hand, appears more evenly dis- 
tributed. All imaging modes in which the nebulosity is 
detected show a dark band centered on the star with a 
position angle of 10-20°. 

The observed nebulosity clearly has a complex 3- 
dimensional structure that cannot be modeled in detail 
by an axi-symmetric code. The model is intended to fit 
an average surface brightness. However, we argue that 
the mid-infrared emission must be optically thin, and 
that the surface brightness therefore is a function of col- 
umn density and UV field only. This allows the model 
to be used to fit an average PAH abundance in the sur- 
rounding cloud. The quantum-heated PAHs are included 
in the model setup as described in Sect. 14.11 The sur- 
rounding material is assumed to have a density profile 




Fig. 4. — Goodness of fit (see text) as a function of PAH abun- 
dance in percent of the the total dust mass and cloud density. The 
contours have been normalized, since the substructure in the nebu- 
losity renders an absolute x 2 value meaningless. The contour levels 
are plotted in steps of 15%. The lighter contours indicate a value 
of 1 and the darkest 3. The constraints used are the IRAC 4 and 
MIPS 1 bands together with the IUE spectrum. The fit is slightly 
degenerate in the sense that a higher density may be coupled with 
a lower PAH abundance, but only as long as the cloud does not 
approach unity optical depth for UV photons within the nebula. 
This happens for densities of > 2000 cm -3 . 

decreasing as Rr 1 , except within 15,000 AU of VV Ser, 
where the cloud density is constant and lowered by a fac- 
tor of 10 relative to the density at 15,000 AU to fit the 
fact that the nebulosity peaks at 15,000 AU. In addition 
to the spherical cavity, the morphology of the nebulosity 
suggests the presence of a conical cavity perpendicular 
to the disk, reminiscent of an outflow cavity. The en- 
velope is thus modeled by the following expression for 
i?>15,000AU: 



where: 



C(0,R) 



P = C{9,R) x (R/Rout)- 



Pout for 6 > (1 - R/Rout) * tt/3.5 
,o out /4 for < (1 - R/R out ) * tt/3.5 



(7) 



(8) 



The parameters of the cavity were not varied in the 
fit. The envelope emission can then be modeled using 
four parameters: the total dust density at 15,000 AU, 
the PAH and VSG abundances and the total extinction 
toward the central star. The extinction or total column 
density of the cloud is considered fixed. This leaves three 
parameters to be varied to fit the surface brightness pro- 
files at 8.0, 24 and 70 /im, where the nebulosity is de- 
tected with sufficient signal-to-noise. 

To find the best-fitting cloud density and PAH abun- 
dance, we calculated a grid of models, comparing each 
model to the imaging data. The model images were 
convolved with theoretical Spitzer point spread func- 
tions, calculated using the STINYTIM package ver. 1.3 
for pre-launch parameters (Krist 2002, SSC web site). 
The goodness-of-fit estimate is the sum of the squared 
differences in pixel values between the observed and 
model intensities in a cross section along the polar axis: 

S[/ bs(2/) ~ ^modcl(j/)] 2 - 

We find that by using only the IRAC 4 band for the 
goodness-of-fit estimate, the model becomes degenerate 
in the sense that a lower cloud density can be compen- 
sated by a higher PAH abundance. This is simply the 
optically thin limit. However, since the relevant optical 
depth is measured in the UV region of the spectrum, in- 
creasing the density too much will cause the nebulosity to 
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Fig. 5. — Cross sections through VV Ser in the Spitzer pass bands oriented perpendicular (left column; positive direction to the west) 
and parallel (right column; positive direction to the south) to the disk plane. The solid black curves are the observations, while the dashed, 
green curves are for the model with the VSGs removed, and the dotted blue curves are for models with both PAHs and VSGs removed. 
The solid, red curves are those of the best-fit model, including both PAHs and VSGs. The model including VSGs is only substantially 
different from the PAH only model outside the disk shadow and at wavelengths longer than ~15^m. At 70 lira, thermal grains dominate 
the emission within 20,000 AU of the central star. 



become optically thick to UV photons at relatively small 
distances to the central star. In the model grid, we find 
that emission above 50 fim, and therefore also the MIPS2 
70 fim emission, is dominated by the thermal grains, 
rather than emission from the VSGs. This provides a 
convenient independent measure of the density of the 
thermal grains. Therefore, by comparing the 8, 24 and 
70 fim bands, the PAH and VSG abundances can be con- 
strained to 5 ±2% and 6 ±3% of the total dust mass. The 
outer cloud density is found to be p out = 500±200cm~ 3 . 
The result is illustrated by the goodness-of-fit surfaces in 
Fig. The best-fitting PAH abundance corresponds to 
2.5 ± 0.8 x IO -5 H -1 or 12 ± 4% of the C atoms be- 
ing bound in PAHs assuming a gas-to- dust ratio of 100 
and a cosmic C/H ratio of 2 x 10~ 4 ( Hol wegen 120011 
lAllende Prieto et alJl2002]) . This is consistent with other 
estimates of the abundance of P AHs in molecular clouds 
and photo-dissocia tion regions l|Boulanger et alJ Il998t 
iTielens et all!9 99). Note that the cloud density matches 
well th at (~ 700 cm" 3 ) determined using the extinction 
map of iCambresvl l)1999T) . 

Fig. shows model cross sections of the best-fitting 
model compared to the Spitzer images parallel and per- 
pendicular to the disk plane of VV Ser. Model cross 



sections for envelopes with no VSGs and with no quan- 
tum heated grains are also shown. This plot illustrates 
the constraints used to fit the envelope. The addition of 
the PAH component is clearly necessary to fit the 8 fim 
image, while the VSGs are required to fit the 24 /jm im- 
age. The 70 /im image shows significant extended emis- 
sion due to thermal grains. The model somewhat undcr- 
predicts the surface brightness in the western part of the 
envelope at this wavelength, although that could possibly 
be corrected by changing the VSG opacity. 

Fig. shows the Spitzer images from 4.5 to 70 fim, 
along with the model fits at 8, 24 and 70 fim. The IRAC 
and MIPS images are seen to be reasonably well fitted in 
terms of the shadow and the south-western emission blob. 
However, the axisymmetric nature of the model makes 
the fit worse regarding other details in the morphology 
of the envelope emission. 

5.2. Disk shadow 

A wedge-like dark band in the mid-infrared nebulosity 
appears naturally in the model and matches well that ob- 
served in the IRAC and MIPS images. This is due to the 
disk blocking UV photons from propagating through the 
disk plane causing a shadow to be cast across the PAH 
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Fig. 6.— Left panels: IRAC4 (8 /an), MIPS1 (24 /an) and MIPS2 (70 /an) images of the VV Ser region. The IRAC4 image is dominated 
by PAH emission and background stars, while the MIPS images are dominated by emission due to a population of VSGs. The nearby red 
object SSTc2d J182850. 21+000949.4 is not related to VV Ser and is possibly a chance projection of a T Tauri star outside the main cloud 
or a background AGB star, right panels: Model contours in the IRAC4, MIPS1 and MIPS2 photometric bands. The contours correspond 
to those of the left side panels. The IRAC4 model contours are overlaid on the SOFI-NTT J-band image to show that there are no visible 
signatures of the mid-infrared nebulosity in the near-infrared. North is up and East is to the left. 



nebulosity where the small grains do not get quantum- 
heated because of the absence of UV photons. While the 
nebulosity and superposed shadow are expected to be 
similar in morphology to those seen in reflection nebulae, 
this is a new type of shadow governed by an entirely dif- 
ferent physical process from single scattering. The open- 
ing angle of the shadow is a measure of the maximum 
height -Htjv /R of the t\jv = 1 surface of the disk. In the 
case of VV Ser, the shadow is wide, corresponding to a 
H\jv/R\jv ~ 0.3. For a self-shadowed disk like VV Ser, 
the maximum H\jy / i?uv is located at the puffed- up inner 
rim. The height at which the inner rim becomes optically 
thick to UV photons is therefore ~ 0.25 AU at the inner 
rim i?i n = 0.85 AU. The inner rim scale height achiev- 
ing this is about 50% larger than that expected from a 
disk in hydrostatic equilibrium. A scale height smaller 
than that of an equilibrium disk would be easier to ex- 
plain, for instance as being due to dust settling. A larger 
scale height is more difficult to understand. Note, how- 



ever, that the height at which the disk becomes optically 
thick to UV photons is degenerate in the scale height 
and the total disk mass, with a larger disk mass compen- 
sating for a smaller scale height. A significant change 
in the disk mass of 1-2 orders of magnitudes is required 
to achieve this effect. Alternatively, since the shadow is 
only probing the vertical density profile at large heights, 
a relatively small change in the wings of the profile at 
large heights while preserving a hydrostatic (Gaussian) 
profile in the disk mid-plane may explain the observa- 
tions. One possibility for such a departure is that small 
grains are lifted to la rge heights by the flow of a disk 
wind ijShu et alJl!994(l . Detailed calculations will be re- 
quired to determine whether this process can produce a 
significant dust column density along the line of sight. 

Also evident is that the shadow is not straight in the 
Spitzer images, i.e. does not follow a straight line through 
the origin, but curves away from this line. For a disk 
shadow, this effect occurs when the disk inclination devi- 
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ates significantly from edge-on and allows the inclination 
of VV Ser to be constrained to 70 — 75° (see also the dis- 
cussion m iPontoppidan fc Dullemondll2005|) . The model 
also reproduces the ratio of surface brightness inside the 
shadow, as compared to that out of the shadow, which 
is an interesting test of the density structure of the disk 
model. 

While the shadow appears in the images because 
quantum-heated grains are excited by high energy pho- 
tons only, the lack of UV photons in the disk plane is ex- 
pected to significantly affect the chemistry of the molec- 
ular cloud as well. In other words, a different chemistry 
is expected to b e seen inside the shado w as compared to 
outside [see also lSemenov et all l)2005|) ]. 

5.3. Near-infrared images 

The near-infrared images from the NTT-SOFI are use- 
ful for obtaining an upper limit to the surface bright- 
ness of any scattered light at these wavelengths. It is 
important that this is consistent with the cloud density 
derived from the mid-infrared nebulosity. The observed 
upper limit on the surface brightness is ^22 and 20.5 
mag/sq. arcsec for the J and Ks bands, respectively. 
The model yields a surface brightness due to scattering 
of J = 23.5, H = 23.1, and K — 22.9 mag/sq. arcsec 
at 20" from the star if the quantum-heating is turned 
off, i.e. these values refer to pure scattering. However, 
including the quantum-heated PAHs increases both the 
J- and iVband surface brightnesses by 2 magnitudes. 
This is due to the near-infrare d continuum opacity intro- 
duced to the PAH opacity in lLi fc Draind (J2001) to fit 
observations of reflection nebulae in massive star-forming 
regions. This means that slightly deeper imaging is ex- 
pected to show nebulosity in the near-infrared. In this 
context, it is also interesting to note that the increase of 
2-3 orders of magnit ude in the near-in frared PAH opac- 
ity compa red to thelLi fc Prairlel l)2001|) opacity as sug- 
gested by iMattioda et al.l l)2005|) would result in easily 
detectable extended emission in the SOFI images. This 
this is not seen, suggesting that the PAH opacity curve 
in the VV Ser nebulosity is similar that of iLi fc Draind 
l)2001t) and that the PAHs are therefore likely mostly neu- 
tral or negatively charged. Alternatively, significantly 
larger PAHs than those used here may also reduce the 
near-infrared extended emission. 

5.4. UV spectrum 

A high quality ultra-violet spectrum of VV Ser ob- 
tained with the International Ultraviolet Explorer (IUE) 
exists (see also Paper I). This provides an important di- 
rect constraint on the UV power injected into the system. 
The UV radiation field in the envelope and the resulting 
PAH excitation can therefore be calculated with some 
confidence. Care has been taken in fitting the SED to 
the UV spectrum. It is interesting to note that the IUE 
spectrum shows shallow absorption from the well-known, 
but unidentified, 2175 A band alway s seen in the g alac- 
tic interstellar extinction curve (see Drainc (2003) and 
references therein). Our model of VV Ser produces a 
2175Afeature 2 times deeper than that observed, indi- 
cating that this band is too strong in the assumed PAH 
opacity. 

Large PAH molecules have previously been iden- 
tified as a possible carrier of the 2175 A band 



(jWeingartner fc Draind 1200 1|) (see also Fig. |5J). How- 
ever, these studies fixed the PAH abundance to fit the 
observed 2175 A band. Since the emission nebula sur- 
rounding VV Ser provides an independent measure of 
the PAH abundance that produces a 2175 A band that 
is even too strong, our results seem to support that the 
PAHs likely contribute significantly to this feature. No 
good fit can be obtained to the IUE spectrum unless the 
2175 A feature is at least dominated by the PAH opac- 
ity. For example, if the 2175 A feature is dominated by 
another carrier, the PAH abundance would be lowered 
significantly (by at least a factor of 4). To compensate, 
a fit to the 8 (im nebulosity would require the density to 
be increased beyond the ~ 1000 cm" 3 where the t\jv ~ 1 
surface passes inside the observed quantum heated neb- 
ulosity. It should be noted that since the actual UV 
opacity of PAHs is not well known, significant uncertain- 
ties remain. In particular, if the actual PAH opacity does 
not ha ve a 2175 A feature, contrary to the ILi fc Draind 
(2001) opacity, then no strong constraints on the PAH 
abundance can be derived from the UV spectrum. 

6. DISCUSSION 

How common are extended quantum-heated nebulosi- 
ties associated with single (low- or intermediate-mass) 
young stars and what do they tell us? In this study, the 
mid-infrared nebulosity around VV Ser is interpreted as 
belonging to the general Serpens molecular cloud, i.e. 
the observed PAH/VSG emission is consistent with a 
constant density cloud. The carrier of the far-infrared 
emission, on the other hand, is not distributed evenly. 
As such the nebulosity is a direct measure of the UV 
radiation field around VV Ser; the nebulosity exists be- 
cause the density (and therefore PAH density) is not 
so low as to not have sufficient PAH molecules, yet the 
density is not so high as to block all UV photons from 
penetrating to the observed edge of the nebulosity at 
30,000 AU. At the same time, there is a strict upper 
limit on the density of small dust grains surrounding VV 
Ser from the lack of reflection nebulosity in the near- 
infrared. The lack of extended emission in the near- 
infrared also puts an upper limit of the near-infrared 
opaci ty of the PAHs, rou ghly consistent with the opac- 
ity of ILi fc Draind l)2001|) . Deeper optical/near- infrared 
imaging may detect the faint reflection nebulosity pre- 
dicted to be there. In the case of VV Ser, this translates 
to a gas density of ~ 500 cm -3 , assuming a dust-to-gas 
ratio of 100 - a value that cor responds w e ll to t hat de- 
rived by the extinction map of iCambresvl l)1999|) . If the 
density is increased by a factor of a few, the PAH neb- 
ulosity around VV Ser becomes significantly more com- 
pact, though still visible, because of extinction of UV 
photons by the large grain population. It is therefore 
expected that Spitzer will detect more of such quantum- 
heated nebulae associated with stars that happen to be 
embedded in the more diffuse parts of molecular clouds. 

Such nebulae may be important for several reasons. 
First, a well-defined stellar photosphere gives the possi- 
bility for accurately calculating the input radiation field; 
the input spectrum can be determined by simply obtain- 
ing a UV-optical spectrum of the star in question. In 
the absence of a star one must resort to using a much 
more poorly defined interstellar radiation field. Using a 
radiative transfer model coupled with a spectrum of the 
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nebulosity, it is possible to directly measure the emissiv- 
ity of the PAH molecules in the cloud for comparison with 
laboratory data. For instance, in the case of VV Ser, it 
can be concluded that the near-infrared PA H opacity is 
low, in accordance with the opacity model of lLi fc Draind 
(2001). Furthermore, by mapping differences in PAH 
features in the radial direction (through a diminishing 
radiation field), it is possible to map the chemistry, hy- 
drogenation and ionization state of the PAHs and couple 
those directly to the radiation field. Finally, the far- 
infrared spectral slope of the nebulosity will constrain 
the composition and size distribution of the VSG com- 
ponent. Such sources may therefore be very important 
for deciphering the role of PAHs and VSGs in molecular 
clouds. 

The fact that in our self-shadowed disk model most 
of the observed disk shadow is created by the innermost 
parts of the disk suggests that the mid-infrared nebu- 
losity may be variable on timescales of 1-10 years (at 
least the dynamical timescale of the puffed- up inner rim) . 
This possible variability should be seen as variations in 
the shadow opening angle. A long term variability of 
VV Ser has indeed been observed in the optical, in addi- 
tion to the UX Orionis events (see Paper I), which may 
suggest that this is reflected also in the shadow morphol- 
ogy. However, it is highly unlikely that the shadow will 
experience strong changes, such as disappearing entirely. 

7. CONCLUSIONS 

Spitzer imaging of the Herbig Ae/UX Orionis star VV 
Ser has shown it to be surrounded by a very extended, 
but bright nebulosity. The nebulosity is detected from 
IRAC band 3 at 5.6 (im to MIPS 2 at 70 /mi. Using the 
disk model for VV Ser developed in Paper I as a starting 
point, we model the extended mid-infrared nebulosity 
surrounding the star by adding small transiently heated 
grains to the surrounding material. 

• We show that the nebulosity is most likely due 
to quantum-heated small grains consisting of both 
PAHs and somewhat larger carbonaceous/silicate 
grains of perhaps ~ 500 atoms. The best-fitting 
abundances are 4-6% and 5-7% by mass for PAHs 
and VSGs, respectively, although these values de- 
pend somewhat o n the opacities used . Assuming 
the PAH opacity of lLi fc Drainel l)2001fl is appropri- 
ate, we find that the 2175 A feature must be mostly 
due to the graphite-like structures within the PAH 
molecules. 

• While the PAHs appear to be evenly distributed 
at distances of 15,000-50,000 AU from the central 
star, the VSG component is more structured with 
most of this material concentrated in a clump ~ 
50" = 12, 500 AU to the south-east of VV Ser. 

• The nebulosity is observed to be bisected by a 
wedge-shaped dark band extending across the en- 
tire nebulosity with a position angle of ~ 10 — 20° 
(east of north) . We interpret this as a large shadow 
cast by the small ~50 AU central disk by blocking 
UV photons from reaching the surrounding cloud. 
This is fully consistent with the disk model devel- 
oped in Paper I. The opening angle of the shadow 



shows that the opening angle of the t\jv = 1 sur- 
face of the disk is H/R ~ 0.3. The inferred disk 
inclination of 72 ± 5° and position angle of (13 ± 5° 
east of north) are important constraints for inter- 
preting the model presented in Paper I, in partic- 
ul ar in relation to the near-infrared interferometry 
of Eisner et all pool . 

• We propose that the VV Ser nebulosity and other 
isolated quantum heated regions excited by a sin- 
gle star with a well-defined UV spectrum may be 
used to strongly constrain the properties of PAHs 
and other types of small grains in molecular clouds. 
For instance, the relation of the ionization state of 
PAHs to the UV field can be accurately mapped 
with a mid-infrared spectral map of the nebula. 

• Additionally, the shape of the mid to far-infrared 
spectrum of the nebula will constrain the size dis- 
tribution of VSGs. In principle, such observa- 
tions coupled with simple radiative transfer calcu- 
lations can serve as a calibration for studies of more 
complex regions containing quantum heated grains, 
such as extra-galactic sources, massive star-forming 
regions and circumstellar disks. 

• Disk shadows cast on quantum heated nebulosities 
are useful for measuring the 7\jv — 1 surface of cir- 
cumstellar disks and therefore provide constraints 
on the grain size distribution in the upper layers of 
disks. We have also shown that due to the large size 
of disk shadows, they are very convenient for de- 
termining inclination and position angles for disks 
that are otherwise impossible to spatially resolve 
through direct imaging. Spitzer-IRAC and MIPS 
imaging surveys are excellent tools for identifying 
quantum heated nebulosities around single stars, 
and Spitzer-IRS is well suited for follow-up spec- 
troscopy. 
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